Divergence of the orbital nuclear magnetic relaxation rate in metals 
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We analyze the nuclear magnetic relaxation rate (l/Ti)or6 due to the coupling of nuclear spin to 
the orbital moment of itinerant electrons in metals. In the clean non-interacting case, contributions 
from large-distance current fluctuations add up to cause a divergence of (l/Ti) H>- When impurity 
scattering is present, the elastic mean free time r cuts off the divergence, and the magnitude of 
the effect at low temperatures is controlled by the parameter ln(/xr), where p, is the chemical 
potential. The spin-dipolar hyperfine coupling, while has the same spatial variation 1/r 3 as the 
orbital hyperfine coupling, does not produce a divergence in the nuclear magnetic relaxation rate. 

PACS numbers: 76.20.+q,76.30.Pk,74.25.Nf 
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I. INTRODUCTION 

The nuclear magnetic resonance (NMR) provides pow- 
erful experimental tool in material science. One promi- 
nent example in the condensed matter is that measure- 
ments and analysis of the Knight shift and nuclear mag- 
netic relaxation rate (also called nuclear spin-lattice re- 
laxation rate) proved to be a decisive tesfci of the validity 
of the BCS theory of superconductivity. In general, the 
NMR signal carries information on the interactions which 
couple the nuclear magnetic moments to their environ- 
ment. In such a situation the theoretical understanding 
of various possible contributions to the total measured 
quantity is very important. 

For metals the most important hyperfine coupling is 
the one with the itinerant electrons^. The magnetic 
hyperfine interaction of each nucleus with the electrons 
may be written as H-hf — ~ 7«I ■ h, where 7„ is the nu- 
clear gyromagnetic ratio and I is the spin of the nucleus 
located at the position R. The units where h = ks = 1 
are used throughout except in Sec. VI. The effective hy- 
perfine magnetic field h has three contributions: 



h(R) = 7e / d 3 r 



8ir _., . s — 3p(s • p) 1 
—sS(p) + ^— --- 

6 p 6 p A 



(1) 



where s and 1 are the spin and the orbital moment of the 
conduction electron located at r, and we use the nota- 
tions p = r R, p = \p\ and p = p/p. In Eq. ([1]) the 
integration is over the sample volume V and 7 e = e/(mc) 
is the electron gyromagnetic ratio, c is the speed of light, 
— e is the electron charge (with e > 0) and m is the free 
electron mass. The first term in Eq. |T]) originates from 
the Fermi-contact hyperfine interaction, the second term 
is due to the spin-dipole hyperfine interaction, and the 
third term is due to the orbital hyperfine interaction. 

In this paper we would like to draw attention to an 
interesting property of the nuclear magnetic relaxation 
rate (l/Ti) or f, arising as a result of orbital hyperfine in- 
teraction: For a perfectly clean metallic system of infi- 
nite spatial dimensions, and in the absence of an external 
magnetic field, (l/Ti) or f, is infinitely large. This fact has 
been already reported by Lee and Nagaosa^. Here we 



provide a thorough discussion of the situation using the 
Green's functions method. The nature of the effect is re- 
lated to the properties of long-range static fluctuations 
of orbital fields and currents. Mathematically, the diver- 
gence of (l/Tx) or & appears upon integrating the relevant 
electronic two-particle correlation function over the mo- 
mentum transfer q: the integral diverges at q = 0. The 
divergence is thus related to the behavior of the electrons 
that are far away from the nucleus. In this paper the in- 
fluence of the finite sample size is included via a cutoff 
of the wave vector integrals, and possible surface effects 
are not considered explicitly. 

The divergence of (l/Ti) or b means that, in principle, in 
a sufficiently clean sample the orbital mechanism for the 
nuclear magnetization relaxation is very efficient and the 
total T\ can be very small. In real material some impuri- 
ties are always present. We show that electron scattering 
off impurities removes the divergence. The magnitude of 
(l/Ti) or b at low temperatures is controlled by the pa- 
rameter In(/Lt7"), where p is the chemical potential and r 
is the mean free time for itinerant electrons. This means 
a logarithmic dependence on the impurity concentration. 
Our numerical estimates for Li and SrzRuOi show that 
even for values of r characteristic of the cleanest samples 
the long-range part of (l/Ti) or 6 is not dominant. But 
it can be big enough, we believe, to be experimentally 
determined. 

The rest of the paper is organized as follows. In Section 
II we present the formula for the relaxation rate which 
is used for the calculation. In Section III we evaluate 
{l/Ti) or b in a simple model in which the normal metal is 
described by the free electron gas. We study both clean 
and impure cases, with the details regarding the vertex 
corrections relegated to the Appendix. In Section IV we 
show that, unlike the orbital hyperfine interaction, the 
long-range contribution to \jT\ due to the spin-dipole 
hyperfine interaction is finite. In Section V we demon- 
strate that the singularity of (l/Ti) or f, in a clean metal 
is in fact quite general. First, we show that placing the 
free electron gas into a periodic potential (the simplest 
model for a crystalline solid) does not remove the singu- 
larity. Next, following the idea of Ref. 0, we show how 
metallic systems with an arbitrary electronic dispersion 
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can be analyzed using the connection between (l/Ti) or j, 
and the nonuniform static electric conductivity <r a/ g(q). 
Recognizing the fact that the static electric conductivity 
of the perfectly clean electronic system is infinite perhaps 
makes the divergence of (l/Ti) or f, in such a system less 
surprising and puzzling. Section VI contains a discus- 
sion, where in particular we estimate the magnitude of 
the effect, and the concluding remarks. 

II. THE RELAXATION RATE FORMULA 

The expression for the nuclear magnetic relaxation rate 
that is most suitable for our purposes has the following 
form£&£: 

_L = -2icoth ImXf_(c ,R), (2) 

where u>o = jnH with H being the external magnetic 
field oriented along z axis. In the derivation of Eq. @ 
the hyperfine interaction is treated as a perturbation for 
the electronic Hamiltonian, which can be quite general. 
The main object in Eq. @ is the Fourier transform of the 
retarded correlator of the effective magnetic fields h(R): 

R) = -ie(t)([h + (t,B.),h-(0,IL)]), (3) 

where h± = h x ± ih y . We take the origin of coordinate 
system at the location of the nucleus (i. e. R = 0). 

The correlator K+_(oj) can be computed starting from 
the explicit expression for the orbital hyperfine fields in 
the formalism of the second quantization. We follow the 
standard procedure and find at first the corresponding 
Matsubara correlator K+_(iv n ) and then apply the ana- 
lytical continuation from the imaginary axis to just above 
the real axis: Kf_(oj) = Kf_(iu n — s- uj + i0+). Alter- 
natively, using Maxwell equations one can express the 
real frequency correlator K+_(uj) through the current- 
current correlator and thus through the electrical con- 
ductivity, for which there exist well developed methods of 
calculation including the kinetic equation approach and 
the linear response theory. 

Based on the magnitude of the nuclear magneton, the 
approximation relevant to experiments in typical labo- 
ratory magnetic fields and not too low temperatures is 
uj -> 0. Then Eq. © gives 

-!- = ->£ km —ImKf _ (w , R = 0). (4) 

The approximation loq — > also implies that statistical 
averages have to be calculated for the system in zero 
external magnetic field H. 

The magnetic field to be inserted in Eq. (0| is the ef- 
fective magnetic field given by Eq. dXJ) . We would like to 
comment that there is one approximation which is im- 
plicit in presenting the fluctuating magnetic field at nu- 
cleus in this form. It consist in neglecting the influence 
of the magnetic field on the electrons themselves. This 



feedback effect may be phenomenologically accounted for 
by replacing h on the left hand side of Eq. ([5]) with 
(1 + 47rx)h, where x is the magnetic susceptibility of 
the material. For ordinary metals \\\ ~ 10™ 5 and we 
will not include this contribution in what follows. The 
complete treatment of the feedback effect, which would 
be necessary when \x\ is large, in particular in super- 
conductors, requires taking into consideration dynamical 
fluctuations of the vector potential A of the electromag- 
netic field. This is beyond the scope of the present paper. 

III. CALCULATION OF {l/Ti) orb FOR A SIMPLE 
MODEL OF NORMAL METAL 

In this subsection we consider in detail the nuclear 
magnetic relaxation rate due to orbital hyperfine interac- 
tion, (l/Xi) or b, in the situation where electronic system 
is modeled by the three-dimensional electron gas (mov- 
ing in the positively charged uniform background). The 
electrons do not interact between themselves, but could 
elastically scatter off randomly distributed nonmagnetic 
impurities. 

This model is sufficient to demonstrate that for the re- 
laxation of nuclear magnetization in the absence of an ex- 
ternal magnetic field (l/Ti) or (, diverges in the clean case. 
This divergence is cut off when impurities are present. 
An extension of the argument to more general situations 
is presented in Section V. 

In the second-quantization representation the operator 
of the orbital magnetic field at the origin R = has the 
form: 

h=-7e / AE^W^^W' ( 5 ) 

which can be obtained from the second term in Eq. ((T|) 
by combining the second quantized representation of the 
momentum operator p with the definition of the angu- 
lar momentum. The hats denoting the operators will be 
omitted below in order to simplify notation. Note that 
expression for the magnetic field is consistent with the 
Biot-Savart law. Indeed, Eq. can be obtained if in 
the expression for the magnetic field 8 at R = 0, given by 
the integral h = — (1/c) J d 3 r[j(r) x r]/r 3 , the standard 
second-quantized representation for the orbital electric 
current is inserted and the integration by parts is per- 
formed while neglecting the surface term. 

A. Free electron gas in the plane wave basis 

The %l> operators can be expanded in the basis of plane 
waves, which are solutions of the Schrodinger equation in 
empty space: 
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where V is the volume of the system. When this expan- 
sion is inserted in the definition of h in Eq. ([5]) a straight- 
forward r-integration leads to the following expression 
for the effective magnetic field due to the electron orbital 
motion: 

k,k' V ; a 



The corresponding Matsubara correlator K¥_{iv n ) 
- /,f dTe iT ^(T T h + (T)h^(0)) has the form: 



K + _{iv n ) = 2 y2 ^ 

ki,k 2 



(k! X k 2 f x + (ki X k 2 )g 

|ki-ka| 4 



x 5(ii/„, ki,k 2 ). (8) 
The quantity 5 is related to the electron bubble: 

t-P 

B(iu n , k^ k 1; k 2 , k 2 ) = / dre <T "- ^ (T T 4 (r) 

xck i . 1 (r)4 2a2 (0) Ck , CT2 (0)), (9) 

where /? = 1/T with T being temperature. The angu- 
lar brackets denote the thermodynamic average with the 
given electronic Hamiltonian, and T T is the imaginary 
time ordering operator. For the free electron gas we have 

S(w„,ki,ki,k 2 ,k' 2 ) = -26 k k >5 k k / ki,k 2 ), 

" " (10) 
and the quantity 5* has the following form: 

m 

_ /(^ )-/(&,) 



£ki - £k 2 



(11) 



where Lo m — 2irT(m + 1/2) and v m = 27rTm, m 
i). ■ 1 are the fermionic and bosonic Matsubara fre- 
quencies, respectively. Also, in Eq. (Ti~Tj) Go(k,iuj m ) — 
\j(iw m — £ k ) is the free electron Green's function, with 

£k = £ k — Mj an d /(^O = l/[ ex P( x /^ 1 ) + l] is the Fermi dis- 
tribution function. Note that S depends on the momenta 
k only through the energies e k = k 2 /(2?n). 

To arrive at Eq. © the factors with the vector prod- 
ucts of four k vectors in the integrand of the expression 
for the correlator K+_ were simplified using the delta 
functions of Eq. (TTU)) : 

(ki x k£) + (k 2 x k 2 )_ (ki x k 2 )2 + (ki x k 2 )2 



(ki-ki) 2 (k 2 -k' 2 y 



|ki - k 2 



For the quadratic dispersion of the free electrons the den- 
sity of electronic states per one spin projection N(e) = 
m\/2me/ (2tt 2 ) does not depend on angular variables. 



Converting the k-summations into integrations and sep- 
arating the angular variables we can write: 



2(47r 7e 



dside 2 N{ei)N(e 2 ) 



x K n (si,£ 2 )S(iv n ,ei,£ 2 ), (13) 



where we have introduced the quantity Kq defined as an 
angular integral, which in the present case contains just 
one nontrivial integration over the angle between ki and 
k 2 that can be performed exactly: 



dO kl dQ k2 (ki x k 2 ) 



a In 



4tt 

a - 



4tt 
1 



ki 



a - 1 



- 2 



(14) 
(15) 



Here fei >2 = |k 1)2 | and a = (k\ + k 2 l )/{2k 1 k 2 ) = (ei + 
£ 2 )/ (2y / ei£ 2 ), which in turn implies that a ± 1 = (ki ± 
k 2 f/{2k l k 2 ) = [ y fEl± s /^) 2 /{2 s /E^. 

Now we perform the analytic continuation of S(iv n ) 
given by Eq. (fTTj) on the imaginary axis to the frequencies 
just above the real axis and expand in lo <C /i, keeping 
in mind that we actually need only the imaginary part 
for small lo. Introducing the notation S(u>) = S(iv n — > 
lo + i0 + ) we write: 

TmS(fij) « 7n5(&-6-o;)[u;/'(6)+0(w a )],(16) 



where /'(£) = df(£)/d£. The delta function in Eq. (Tj 
allows us to eliminate one energy integral in Eq. (| 13[) . 
while the derivative of the Fermi distribution places the 
energy close to the Fermi surface. This allows us to use 
for the degenerate electron gas the standard approxima- 
tion of the constant DOS at the Fermi level and switch 
in Eq. (|13|) from the e integrations to the £ integrations 
from — oo to +oo. Additionally, the angular factor Kq 
reduces to 



l?l-?2| 



(17) 



where p, = (4/e)/z « 1.47/x. It is worth emphasizing that 
the angular integral Kq brings in a nonanalytic energy 
dependence. Thus, at small lo we arrive at the following 
integral: 



, TmK2_(w) 2ir 
hm — — 

uj^O LO 



2-k r +oa 

= —[A-K le N{n)f J dfcdfc/' (6) 



x Hti -6 - In- 



/' 



(18) 



[47r 7n7e iV( A1 )] 2 ln(£) , (19) 



16-61 

The final expression for the relaxation time has the form: 
1 _ 2tt 
J\f ~ T w, 

where the right-hand side does not depend on tempera- 
ture. Note that this expression diverges logarithmically 
at to — > 0, which can be traced back to the logarithmic 
behavior of the angular integral Kq in Eqs. (fTS")) and (fTJ 
at £i — > e 2 or ki — > k 2 . 



4 



B. Free electron gas in the spherical harmonics 
expansion 

It is instructive to trace the origin of the divergence us- 
ing the basis of spherical harmonics around R = 0, where 
nuclear spin is located. For this purpose we substitute in 
Eq. ([5]) the following expansion of the plane waves: 

exp(zk • r) = 4tt^ ]T i l 3i{kr)YC m {^)Y lm {v). (20) 
;=0 m=— I 

in terms of the spherical harmonics Yi m and the spherical 
Bessel functions ji(kr)^. An advantage of this approach 
is that we can use in Eq. ([5]) the well known matrix ele- 
ments of the angular momentum operator 1. The radial 
part can be integrated to give: 



(47T) 2 _ f ~ 



V 



kik^cr 



1=1 



4fcifc 2 



(fci + k 2 y 



V(l ± m)(l + m + VjY^MY^fa), 



Zi(z) 



T(l) 



4 r(Z + 3/2) 



-J 1 + 1; 2(1 + !);«], 



where k = k/|k| and F[a, b;c]z\ is the hypergeometric 
function^. Calculation of the Matsubara correlator pro- 
ceeds as in the previous subsection and produces S{iv n ) 
of Eq. (fTT]) and S(us) of Eq. (jTHJ) after the analytic con- 
tinuation, leading to <5(eki — £k 2 ) an d k\ = k 2 &t lj = 0. 
At this stage the angular parts of k-integrations are per- 
formed and then all the summations except one can be 
evaluated. The final result reads: 



1 



2tt 



2 E 

i=i 



16/tt 



(21) 



where finite terms have been omitted. Since we have set 
up explicitly u> = this expression is divergent. In accor- 
dance with what has been found previously in Eq. (|19p . 
the divergence is logarithmic and in this case occurs in 
the sum over orbital quantum number as I — * oo. 



C. The impurity dependence of (l/Ti) or 5 

In this subsection we consider the electron gas inter- 
acting with point-like elastic scatterers. The disorder is 
treated perturbatively using the standard impurity av- 
eraging technique 10 . In the Appendix we show that in 
the isotropic situation considered here the vertex correc- 
tions to the impurity averaged electronic bubble, given 
by Eq. ([9]), vanish exactly. Therefore, the impurity aver- 
aged correlator of the orbital magnetic fields K¥_(iv n ) is 



given by the same expression as in the clean case, Eq. ([8j , 
with unchanged angular part Kq 7 Eq. (| 15(1 but with 



S{iv n , ki, k 2 ) = G(ki ,iui m + iv n )G(k 2 ,iui m ), (22) 



expressed through the impurity averaged Green's func- 
tions: 

G(k,«w m ) = [iuj m - £k + isgn(w m )/(2r)] _1 . (23) 

Here 1/(2t) = nN(fi)ni mp V^ np is the elastic scattering 
time where 7ij mp is the impurity density and Vi mp is the 
Fourier transform of the impurity potential. The Mat- 
subara summation in Eq. (|22|) can be conveniently per- 
formed using the spectral representation for the Green's 
functions. The result is: 



lim 



ImS(cy,6,6) 



to' 



+ OC 



dwf'{uj)A{^,uj)A^ 2 ,uj), (24) 



where A(£,ui) = A(^,u)) = — ImG(k, iu) m — > uj + {0 + )/tt 
is the electronic spectral function in the dirty normal 
metal which is sharply peaked around £ = uo. The 
presence of f'(u>) in Eq. (|2~4"|) then allows us to switch 
to the £ integration with the infinite limits in the ex- 
pression for the correlator K+_(u) as in the clean case 
(see Eqs. (fH?)) and (J5]) ). Now these £ integrations give: 

S±~dS 1 dt a A(ui,t 1 )A(ui,t a )K a (t 1 ,&) = |ln(TAi), which 
leads to the following final result: 

jTf = Y [^InleNip)] 2 In (r/2) . (25) 

We see that 1/r replaces the frequency w, which was 
present in the correlator for the clean case (see Eq. HH)), 
and thus removes the divergence at u> — > 0. 



IV. COMPARISON OF THE ORBITAL 
CONTRIBUTION TO 1/Ti WITH OTHER 
CONTRIBUTIONS 

We extend our analysis to include all the hyperfmc 
magnetic fields of Eq. (TT]) , using the same model for elec- 
trons as in the previous section. In the plane-wave basis 
[see Eq. (|6])] one has for three-dimensional system: 



47T 7 e 

' 3 V 



E 



k,q 



k,q 



q(q • cr a p) <j a p 



V ^ 



k.q 



q x k 

—2 



" C ka Ck -q ! "' 



(26) 

4c k _ q ,^7) 
(28) 



In these equations the summations over all repeated spin 
component indices are assumed. We observe that the 
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operator of the Fermi-contact effective magnetic field, 
h^_ c , and the operator of the spin-dipole effective mag- 
netic field, h s -d, can be combined into a single object, 
for which the q-dependence inside the sums appears only 
through the transverse projector Sij — qiqj/q 2 . This sim- 
plifies the next step of the relaxation rate calculation, 
which is to evaluate the Matsubara correlator K¥_ . We 
obtain: 



X-(-n) = 2 (2?r7 ' 



V 2 

q x k) + (q x k 



E 

k,q 



2 - 



q+q- 



SO„,k,k-q). (29) 



In the clean case the quantity S is given by Eq. ([TT]) . In 
Eq. (j2"9")) the first two terms in the square brackets repre- 
sent the combined contribution of the Fermi-contact and 
the spin-dipole hyperfine interactions. The third term is 
the contribution due to the orbital hyperfine interaction, 
which was considered before. It is the same as in Eq. ([5]) 
but written with the change of the summation variables 
(ki,k 2 ) -> (k,q) where k = (k : +k 2 )/2, q = k x - k 2 . 
In this form the behavior of different contributions as 
functions of q = |q| is clearly displayed. In fact, both 
Fermi-contact and spin-dipole contributions depend on 
q through S function only, while the orbital contribution 
contains the additional factor l/q 2 . 

The divergence in (l/Ti) or (, in the clean system ap- 
pears now as the divergence of the q sum as q — ► 0. 
On the other hand, the spin-dipole hyperfine interac- 
tion produces a finite relaxation rate. This result contra- 
dicts expectations 4 that spin-dipole and orbital hyperfine 
interactions should bring about similar contributions to 
1/Xi because in the real space they both vary with the 
distance from the nucleus as 1/r 3 [see Eq. (JTJ]. 



V. GENERAL CHARACTER OF THE 
DIVERGENCE OF (l/Ti) orb IN CLEAN METALS 

In this section we demonstrate that the divergence of 
I/Ti in a clean metal is quite general and is not an ar- 
tifact of the free electron model we used in previous sec- 
tions. We now concentrate on the orbital hyperfine in- 
teraction that is responsible for the effect. First we con- 
sider the electron gas in an arbitrary periodic potential, 
and show that the divergence does not disappear. Then 
we relate the correlator of orbital magnetic fields to the 
current-current correlator, and therefore to the electri- 
cal conductivity. This formulation allows one to treat an 
arbitrary electronic dispersion quite generally. The di- 
vergence is shown to be a consequence of the behavior 
of the static nonuniform conductivity in the long wave- 
length limit. 



A. The electron gas in a periodic potential 

Now add to the Hamiltonian of the clean Fermi gas 
a periodic potential, and consider the effect it has on 
the correlator of orbital magnetic fields in the definition 
of I/Ti in Eq. ([2]). The solutions to the corresponding 
Schrodinger equation are the Bloch functions which can 
be used as the basis for the second-quantized description: 



^EE< 

v n k 



/kr 



w„k(r)c„ : 



k<T J 



(30) 



where n enumerates the electron bands and the wave vec- 
tor k is now limited to the first Brillouin zone. Func- 
tions M n k(r) are periodic in r and can be expanded in 
the Fourier series: 



ink(r) = ^[/, iG (k)e iC 



(31) 



where G is a reciprocal lattice vector. The operator of 
the effective magnetic field at the origin, R = 0, due to 
the electron orbital motion has the form: 

Airi \ - \ - (k + G) x (k' + G') 
- 7e ~ (k + G-k'-G') 2 

kk' GG' v ' 



X 



J2 U nG (k)*U n , G ,(k') 4k CT <VkV- (32) 



This expression is to be compared with Eq. ([7]). The 
Green's functions of the electrons in a periodic potential 
are diagonal in k as well as in the spin and band indices. 
Using these properties we calculate the correlator of the 
orbital magnetic fields at the location of the nucleus. The 
expression for the relaxation time in the limit uj — > can 
be written as follows: 



1 (47T7»7 e ) ^ ^ 

Iff = 27T — <H £ ™i*l - £n 2 ,k 2 ) 

1 " kik 2 ni"2 



X 



/'(£n 2 ,k a ) ]>> 4 K 1+ K 2 _ 

{G} 



(33) 



where {G} = {G u G 2 , G 3 , G 4 } and 



W = U. 



nl,Gi + G : 



.(kijC/^cJki) 



xU* 



Ki = 
K 2 = 



y ri2,G 3 +G 4 (M 2 )U n2 ,G 2 (k 2 ), 

(kg ~k 2 + Gi) x (k 2 + G 2 ) 
(ki-ka + d) 2 
(k 2 - ki + G 3 ) x (ki + G 4 ) 
(k 2 - ki + G 3 ) 2 



(34) 
(35) 

(36) 



Eq. (|33[) is an exact expression. We show that it con- 
tains an infinite term. To separate the divergent contri- 
bution we restrict the sums in Eq. (|33[) by the conditions 
Gi = G 3 = and n\ = n 2 , and change the summa- 
tion variables ki = k + q/2, k 2 = k — q/2. Then we 
perform the small q- expansion of the arguments of the 
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delta function and /' in Eq. ((33)) , as well as of the U 4 of 
Eq. (f3~4")> . Accounting only for the leading term in these 
expansions Eq. (|33|) is reduced to: 



TiT 



V 2 



k,q n 



[q x v n (k)] + [-q x v„(k)]_ 



(37) 



where we have used the following expression^ for velocity 
of the band electrons: 

V «W - ^E( k + G )I^G(k)| 2 = (38) 

G 

To estimate the q- sum in Eq. (|3"7f we change it, for each 
band n, to the integral over the sphere of an appropri- 
ate volume. The integration over q = |q| extends from 
q m m = to q max — 2q Fjn , where q F , n is of the order of 
the average separation between two points on the Fermi 
surface in the n band. The final result has the form: 



k.n 



K(k)] 2 + ^Vn(kf 

|v„(k)| 



(39) 



where D n = hi(2qF,n/qmin) is the divergent factor. As 
long as the electrons do not interact with each other 
q-min = at any temperature. The reason is that the 
delta function, which appears in Eqs. (f3"3"| and ([57)1 and 
ultimately brings in the divergence, is in fact the elec- 
tronic spectral function of the noninteracting system. 
Thus, the same logarithmic divergence as found in Sec. 
Ill for free electrons appears here. For finite samples of 
the linear size L we have q m in = 1/L. When impurities 
arc present, from the result of Sec. Ill C it is expected 
that qmin — 1/(vft), which is in general different for 
different bands. In the next subsection we confirm this 
form of q m in for one band with an arbitrary electronic 
dispersion. 

To conclude this subsection we note that for band elec- 
trons we found the logarithmically divergent contribution 
to (l/Ti) or b and the factor in front of it (see Eq. l[3"3"))). 
which depends only on the averages of the Fermi veloci- 
ties and does not involve the position of the nuclear spin 
in the unit cell. 



B. Expressing K+_(uj,q) in terms of the 
current— current correlator 



We now use another method to demonstrate the singu- 
lar behavior of the nuclear magnetic relaxation rate due 
to orbital hyperfine interaction. It also allows us to ob- 
tain a useful expression for (l/Ti) or {, for an anisotropic 
metal. The idea, first used by Lee and Nagaosa^, is 



to express the correlator of effective orbital magnetic 
fields K+_(uj 7 II = 0) = ^2qK^_(u>,q), which appears 
in Eq. (J3)), through the current-current correlator, and 
consequently through the nonuniform electrical dc con- 
ductivity. 

The fluctuating magnetic fields are calculated using 
the Maxwell equation, assuming that the electric cur- 
rents are given: V x h(u>, r) = (47rj(o;, r) — iwE(w, r)) /c. 
In metals it is a good approximation to neglect the sec- 
ond term containing the electric field E. Then, applying 
the Fourier transform to find the magnetic field we get 
h(w,q) = (47ri/c)[q x j(cj, q)]/<? 2 , and the required re- 
tarded correlator of the orbital magnetic fields is given 
by the following expression: 



(A-K\ 2 q m q k r 



xmi^xkj T" Cymi^ykj 



ie z ipei mi e p kj]llij(uJ,q), (40) 



where Hij(u>,q) is the retarded current-current correla- 
tor, eijk is the totally antisymmetric tensor and summa- 
tion over the repeated indices, which run through x, y 
and z, is assumed. In the end we are interested in the 
situation when the constant external magnetic field is ab- 
sent. Therefore we take the 11^ tensor to be symmetric, 
namely Ily = Hji. In this case the third term in the 
square brackets in Eq. (|4T)]) vanishes. 

In general ELy contains both paramagnetic and dia- 
magnetic terms. However, for the calculation of 1/Ti we 
need only Imlly (see Eq. fl2J)), and the diamagnetic part, 
which is real, drops out. The retarded current-current 
correlator can be related to the electrical conductivity 
tensor—. Using ImlXy (w, q) = — Re[wtTy (ui, q)] we ob- 
tain from Eq. (|4"0")) the following expression: 



lxnKf_(uj,q) _ /4tt\2 1 
-—II U V c / q 4 

,2, 



lim 



^e[(ql+ql)a z 



+ q z {<?xx + (Tyy) - 2q x q z a xz - 2q y q z a yz \ , (41) 

where arguments of the conductivity tensor <7ij(uj = 0, q) 
are not shown explicitly. The expression (|4"T)) has to be 
inserted in Eq. ((4]). 

In the simple case of an isotropic system the conduc- 
tivity tensor is diagonal in the coordinate frame with 
the z axis parallel to q — (q,9,<p), namely Re cry (g) = 
di&g[a±(q), cr±(q), au(q)]. This tensor should be trans- 
formed to the coordinate frame defined by the initial ori- 
entation of the nuclear magnetic moment (i.e i||H) using 
an appropriate rotation of coordinate frame and then in- 
serted into Eq. (|4T)) and Eq. (j4)) . The final expression for 
nuclear magnetic relaxation rate contains only the per- 
pendicular conductivity a±(q): 



1 



dqa±(q). 



(42) 



It is known from the theory of the anomalous skin effect^ 
that in a clean metal the nonuniform conductivity <jj_ (q) 
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behaves as 1/q at small wave vectors q. This leads to 
the logarithmic divergence for the nuclear magnetic re- 
laxation rate (l/Ti)or&- 

In order to clarify this statement and obtain a useful 
formula for (l/Ti) or {, for electrons with an arbitrary dis- 
persion in the presence of impurities, the q-dependent 
electric conductivity tensor is now evaluated using the 
kinetic equation approach. Electric current is given by 



j(w,q) = 2e^v k n(u;,k,q), 



(43) 



where v = Vk = <9£k/9k is the electron group velocity, 
and the factor 2 is due to the spin. The distribution 
function n(oj, k, q) = /(^k) + 6n(oj, k, q) is a sum of the 
equilibrium Fermi distribution / and a correction Sn pro- 
portional to the small perturbing electric field E(q) = E. 
Within the relaxation time approximation for the impu- 
rity collision integral we have: 

[1/t - i(u - q- v)] Sn = -e(df/d£)v ■ E, (44) 

leading to 



5n(ui, k, q) 



eT(df/d£)v ■ E 



1 — ir(u — q • v) 
The final result for the conductivity reads: 



(45) 



fjy(w,q) = 2e 5 



T (1 + ITUj) V l V° 

1 + t 2 (uj — q • v) z 



47T 



(46) 



This form of q) is valid for an arbitrary electronic 

band £k, which enters through the angle dependent den- 
sity of electronic states N(£ , f2k). Using the form of the 
conductivity tensor form Eq. (|46[) in Eq. (|41[) we arrive 
at the following expression for 1/Ti due to the orbital 
hyperfine interaction: 



1 



8 /e 7 n\ 2 f°° dq [°° ( df 

2(fe^)(q- v) 



v 2 z + v 2 q 2 



l/r 2 + g 2 (q-v) 2 



(47) 



where the z direction is set by the initial orientation of 
the nuclear polarization. The angular brackets denote 
averaging over both k- and q- angular variables. The 
magnitude of the group velocity v = Vk is fixed at the 
Fermi velocity by the derivative of the Fermi distribu- 
tion. For an anisotropic system the Fermi velocity is an- 
gle dependent. For electrons with quadratic dispersion 
the calculation of 1/Ti based on Eq. (|47|) gives the result 
identical to that given by Eq. (j42|) . 

The q integral should be cut off by 2q F at the up- 
per limit, which is necessary in the present approach due 
to the fact that the kinetic equation does not capture 



the large-q behavior of the nonuniform conductivity cor- 
rectly. On the other hand, the q integral is convergent at 
the lower limit. The clean limit corresponds to 1/r — > 0. 
In this case the right hand side of Eq. (|4T|) contains the 
following 5 function: 

1/a K T , a ^^(q-v)] = -i[q-v]. (48) 
l/r- + (q-v)" : q 

After the integration over 17 q is performed we can see 
explicitly that the remaining q integral diverges as 1/q 
at the lower limit. 



VI. DISCUSSION 

It was found experimentally 12 that the nuclear mag- 
netic relaxation in many metals is dominated by 
the the Fermi-contact hyperfine interaction, for which 
(l/T 1 T) J r_ c = (Ank B /9)[47r 7n7e h 2 N(n)] 2 (\u k (0)\ 2 ) 2 FS 
[in this Section we restore k B and H]. The situations 
where this is not so and the contributions of hyperfine 
orbital and/or spin-dipole interactions are significant are 
worth of detailed analysis. Loosely speaking, the physi- 
cal reason for the large magnitude of the Fermi-contact 
interaction is that the conduction band usually contains 
a large portion of the atomic s-orbitals which give a large 
contribution to the electronic spin density at the nuclear 
site, specified by the overlap factor (\uk(0)\ 2 }' FS . The 
relative importance of the other hyperfine interactions 
increases when the density of states at the Fermi level is 
dominated by the d— and /- bands as in transition metals 
and their compound o 13 i 14 ' 15 i 16 , or when the conduction 
band happens to be almost free of s orbital o 17 i 18 ' 19 . 

We would like to emphasize that both orbital and 
spin-dipole hyperfine contributions to 1/Ti commonly 
discussed in the literature are in fact the "local" contri- 
butions which we did not touch upon in this paper. It 
is usually argued that since the orbital and spin-dipole 
hyperfine interactions decay fast enough with distance 
[they both vary as 1/r 3 , see Eq. (JT])], the nuclear spin 
of a given atom is only affected by the magnetic fields 
generated by the electron orbitals centered at this very 
atom. In contrast, the contribution we analyzed in the 
previous sections comes from large distances. It has to 
be added to the local orbital contribution. 

Below we attempt to estimate the magnitude of the 
relaxation rate due to the long-range part of the orbital 
hyperfine interaction [denoted (1/Ti)^] in some metals, 
and compare it with the contributions due to local part 
of the orbital hyperfine interaction [denoted {1/T\) l °^\ 
and Fermi-contact interaction. We first consider alka- 
line metals. In this case, the estimates based on the 
quadratic dispersion produce the correct order of magni- 
tude for (1/Ti)ir_ c , as compared to experimental values. 
This makes our further estimates of (1/Ti)^ trustwor- 
thy. Next we discuss Sr2Ru04, the compound which is 
likely to be suitable to study the orbital nuclear magnetic 
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material 


T 1 T F-c 


i 

orb 


1 orb 


(1/T!) F _ C 




Li 
Na 
Sr 2 RuOi 


0.019 
0.213 
1.L10" 4 


4.05-10" 4 
1.20-10" 4 
1.3-10" 5 


1.14-10" 3 


8.5-10 -4 


2.2-10 -2 
5.6-10 -4 
0.12 


0.37 
0.015 



TABLE I: Numerical values of various contribution to 
l/(TiT) [in the units (K • sec) -1 ] for some metals obtained 
as described in the text. For Sr2RuOi the data are related to 
17 in the position 02, which is out of the conducting basal 
plane. 



relaxation due to the peculiarities of its electronic struc- 
ture. Additionally, in this case very clean samples with 
a long mean free path are available. The results of our 
estimates are summarized in Table Q] 

A. Alkaline metals 

We consider JVa and Li metals. For both of them the 
Fermi surface lies within the first Brillouin zone and the 
electronic dispersion can be taken e(k) = /7, 2 k 2 / (2m*) for 
all energies up to the Fermi level 2 ^. This means that the 
results of Section III C are directly applicable. Numerical 
values are m*/m = 1.45, q F — 1.3 ■ 10 8 cm -1 for Li, and 
m*/m = 0.98, q F = 9.2 • 10 7 cm -1 for JVa, as quoted 
in Ref. I2H The nuclear magnetic moments are 3.26/ijv 
for 7 Li and 2.22/ijv for 23 Na, where /ijv is the nuclear 
magneton. 

In the case of the Fermi-contact interaction the neces- 
sary overlap factors were obtained using the results of the 
first principle calculations of Kohn 22 and Kjeldaas and 
Konh 23 : 16 2 for Li and 144 2 for Na. Thus, our theoretical 
estimate for Li is {l/TxT) F - c m 0.019(AT • sec) -1 , which 
is to be compared to experimental values 0.023(N-sec) _1 
of Ref. [24| or 0.06(K ■ sec)" 1 of more recent work in 
Ref. In the case of Na our theoretical estimate is 
(l/TiT) F _ c « 0.213(fsr • sec)- 1 , while the experimental 
value from Ref. H is 0.196(N • sec) -1 . 

The expression for the nuclear magnetic relaxation 
due to the long-range part of orbital hyperfine inter- 
action for quadratic dispersion can be obtained from 
Eq. ((23 and Eq. @g) as follows: (1/TiT)^ = 

(2TT/3)[i7Tj n j e h 2 N(^)] 2 (m/m*) 2 ln(2q F l). For our esti- 
mates we use I = 10 -5 cm and obtain that {^-/TiTyZT = 
4.05 • 10 -4 (if • sec) -1 in Li and 1.20 • 10 -4 (if • sec) -1 in 
Na. We present here a useful formula that defines the 
ratio of the two discussed relaxation rates: 

(1/T X ) F _ C 2\m*J (K(0)| 2 ) 2 S - ^ ' 

For Li this ratio is 2.2% while for Na it is 0.056%. Such 
large difference is explained by the ia,ci^ that the wave 
function of Li at the Fermi level contains about 73% of 
the I = 1 spherical harmonic, and only 20% of the spher- 
ically symmetric I = harmonic. 



The magnitude of the local part of the orbital mag- 
netic nuclear relaxation in both metals was obtained 
experimentally by Hecht and Redfield 25 from the mea- 
surements of the Overhauser effect. They found that 
[(l/Ti)^ c b ]/[(l/Ti) F _ c ] ps 0.06 ± 0.03 in Li, while there 
is no local contribution in Na. Thus, the local and long- 
range contributions to l/(Ti) OT -6 appear comparable in 
Li, and both are quite noticeable. 

It appears that metallic Li could be the system of 
choice to look for the long-range orbital contribution to 
the nuclear magnetic relaxation. One can use the char- 
acteristic logarithmic dependence of {1JT\T)^ on the 
impurity concentration. 

It is worth noting that the NMR experiments are 
often done with collections of relatively small metallic 
particles . This is dictated by the experimental require- 
ment that a large volume of the sample should be sub- 
jected to a radio frequency magnetic field which only pen- 
etrates within the skin depth. Such experimental setup 
would limit the magnitude of the long-range orbital ef- 
fect even for the very clean metals. 

B. Sr 2 Ru0 4 

The electronic properties of this material are well 
studied 2 - 7 -, both theoretically and experimentally, in con- 
nection with unconventional superconductivity, which 
appears at temperature of about IK. For the purpose 
of a rough estimate for the long-range orbital magnetic 
relaxation rate we use the model of three two dimensional 
bands, with cylindrical Fermi surfaces and q F ,i = 3.04 • 
10 7 cm -1 , q F:2 = 6.22-10 7 cm" 1 and q F , 3 = 7.53-10 7 cm -1 , 
as quoted in Ref. The lattice parameter along the 
tetragonal symmetry axis is 03 = 1.27 ■ 10 cm. To ob- 
tain a convenient formula from Eq. (|39p we assume that 
v z/( v x + v y) ^ 1 an d the cutoff wave vector q min is the 
same for all bands. Using the density of the electronic 
states for each band in the form N 2 d(a 1 ) = m* / (2irh 2 a?) 
and taking z along the tetragonal symmetry axis we ob- 
tain: 

^" = ^P) 2 E^M2^J). (50) 

T 1 T orb h V C 1 t^l ° 3 

The best samples of Sr-iRuO^ used in the de Haas - 
van Alphen experiments 2 ^ are very clean, with the elas- 
tic mean free path reaching 3 • 10 -5 cm. For the re- 
laxation of the magnetic moment of 17 nucleus, with 
7„ = -3.63 • 10 3 (sec • G)"\ we obtain (1/Tj.T)^ = 
1.3-10~ 5 (-Kf -sec) -1 . This value can be compared with the 
local part of the orbital nuclear relaxation rate. For ex- 
ample, from the band structure calculations Pavarini and 
Mazing found that (1/TiT)^ = 8.5 -10 -4 (if • sec) -1 for 
the so called 02 oxygen, located out of the basal conduct- 
ing planes. The ratio is [(l/Ti)^] / [{l/T x ) l ™ b ] =0.015. 

The 02 position of the oxygen is very favorable for 
studies of the long range orbital mechanism of the nu- 
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clear magnetic relaxation, because the local orbital con- 
tribution is much smaller than for the other oxygen site, 
01. The Fermi contact contribution is also numeri- 
cally quite small for the 02 position, (l/TiT)p _ e = 
1.1 • 10~ 4 (if • sec) -1 as calculated in Ref. 19. It is still 
larger than the long-range orbital contribution, but not 
substantially: [(1/Ti)£j] / [{1/Tx)f- c ] =0.12. 

Comparison of our predictions with experiments is 
complicated by the fact that the experimental values 
of 1/T\T for 02 oxygen 2 ^ are two orders of magnitude 
greater than values given by the band structure cal- 
culations. The experimental results necessarily include 
the effect of electronic correlations, which are strong in 
S12R11O4. Such effects were completely ignored in our 
discussion of (l/Ti)^ in this paper. 



C. Concluding remarks 

We would like to point out that besides the elastic im- 
purity scattering, other scattering mechanisms for elec- 
trons should also affect to (l/Ti) or f,. One expects that 
inelastic scattering, for example by phonons, should lead 
to 1/(T\T) ~ InT, because of the temperature depen- 
dence of the inelastic scattering rate. 

The NMR relaxation mechanism through the orbital 
hyperfine interaction that we have discussed in this pa- 
per is quite general. It is not limited to metals and should 
exist in any system with mobile charge carriers. On the 
basis of Eq. (|19")) one expects that a small effective mass of 
the charge carriers should significantly enhance the effect. 
Small masses are typically found in semiconductors 2 ^. In 
such systems the NMR techniques have also been very 
successful, though the analysis should be modified com- 
pared to the case of metals^. 

Finally, we would like to comment on calculations of 
(l/Ti) or b from the first principles using band structure 
methods. These calculations typically use the expansions 
in a set of basis functions that involve the spherical har- 
monics centered at the nuclear site, and an appropriate 
radial part obtained from the solution of the Schrodinger 
equation inside the Wigner-Seitz sphere or the muffin-tin 
spher o 15 i 16 ' 18 . The series in the orbital quantum number 
I is truncated at some stage of the computation based 
on an empirical estimate of convergence. However, we 
have demonstrated in this paper that this expansion is 
expected to diverge as In I, which is usually difficult to 
capture numerically. The presence of this divergence calls 
for a careful examination of the conditions used to justify 
the truncation of the I expansion. 

In conclusion, in this paper we have discussed the nu- 
clear magnetization relaxation rate due to the orbital 
hyperfine interaction of nuclear spins with itinerant elec- 
trons. For an infinite clean metal (l/Ti) or h would become 
divergent at zero temperature. The reason for the diver- 
gence is that the total fluctuating magnetic field created 
by all the electrons at the site of a given nucleus contains 
a contribution from distant electrons. The divergence 



is removed when there are scattering mechanisms in the 
system. At low temperatures the scattering is predomi- 
nantly elastic. In this case (l/Ti) or b depends logarithmi- 
cally on impurity concentration. Based on our estimates, 
it seems to be feasible to observe this effect e. g. in metal- 
lic Li. 
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APPENDIX A: IMPURITY AVERAGED 
CORRELATOR OF ORBITAL MAGNETIC 
FIELDS 

The operator of the magnetic field created by the or- 
bital motion of the electrons at R = is given in Eq. {7J . 
Here we want to compute the impurity averaged correla- 
tor of such orbital magnetic fields at a given Matsubara 
frequency iv n . This quantity can be written as follows: 

Ik-m c u _ (47r 7e ) 2 x ^ (ki x k;) + (kg x k^)_ 

[K + _{lV n )) imp - y2 2^ (kl _ k /)2 ( k2 _ k > )2 

x (B(iv n ,k 1 ,k' 1 ,k2,k' 2 )) mp: (Al) 

where {k} = ki,ki,k2,k 2 and {...)imp means the aver- 
age over impurity configurations. The electronic bubble 
B(iv ni ki,k' 1; k 2 , k' 2 ) is defined in Eq. ((9]). For the follow- 
ing discussion it is more convenient to use a different set 
of wave vector variables, namely Q a = (k a + k' a )/2 and 
q a = k a - k' a with a = 1, 2. 

The impurity averaging restores the translational in- 
variance of the system and the result for the electronic 
bubble, after the spin summation is performed, can be 
written as 

(Z?(w n ,Qi,qi,Q2,q 2 ))imp = -2<5 qi -q 2 S Ql: Q 2 (iu n , qi). 
The function S has the following general structure^: 

< 5 Qi,Q2(M'n>q) = T G(Qi+, icj m+ )G(Qi_, iwm) 

m 

x [*Qi,Qa + 4 r Qi,Q2(q)* w TO>^n) 

x G(Q 2+ ,iw m+ )G(Q 2 -,KJ m )], (A2) 

where Q a± = Q a ± q/2 for a = 1,2, uj m+ = u m + v n 
and G(k, iui m ) is the impurity averaged Green's func- 
tion defined in Eq. (|23p . The first term in the square 
brackets in Eq. (|A2J> corresponds to the bare bubble con- 
tribution to S. Here the interaction with impurities is 
included through the electronic Green's functions only. 
The second term in the square brackets in Eq. (|A2|1 repre- 
sents the vertex corrections, which are defined in terms of 
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the four point vertex function Tq^c^ (q, iu> m , iv n ). This 
function satisfies the following integral equation 31 : 

1 x - 

r Qi,Q 2 = U q^q 2 + 7/2^ U QuQ 
Q 

x G(Q + ,ico m+ )G(Q-,iuj m )TQ t Q 2 , (A3) 

where the dependence of T and U on q, iuj m and iv n is 
not shown explicitly. 

The function J/q^q,, (q, iuj m , iv n ), which appears in 
Eq. (|A3j) . contains all the information about the effec- 
tive interaction between electrons arising due to impurity 
scattering, and can be very complicated in general. For 
good metals it is usually a good approximation to con- 
sider impurities as point-like scatterers and take them 
into account in the Born approximation. In this case 
the function C/q 1j q 2 (q, iu> m , iv n ) is essentially a constant, 
equal to u = n imv Vf mp . Then Eq. (|A3|) has the following 
simple solution: 

T~ l (q,iu) m , iu n ) = - - i V* G(Q+, iw m+ )G(Q_, iu) m ). 
u V 

Q 

This result contains the sum of the ladder diagrams and 
correctly reproduces the diffusive dynamics of the elec- 
trons at low wave vectors and frequencies. 

We return now to the correlator of the orbital mag- 
netic fields Jf| f _(w„) in Eq. (|A1|) . and split it into the 
two pieces: K+_(iv n ) = K bare (iv n ) + K vert (iv n ) cor- 
responding to the two terms in the square brackets of 



Eq. (|A2[) . The bare bubble contribution has the form: 



K bare {iu n ) = 2 [ ™ (q x Q)+ (q x Q)- 
x-j^Q.q,*^), 



where the function S was defined in the main text in 
Eq. (|2"2")1 , but using the variables ki = Q + q/2, k 2 = Q — 
q/2. This bare bubble contribution has been analyzed in 
Section IIIC. The vertex correction part can be written 
as: 



K vert (iu n ) = 2^T^(qxL(q)) + (qxL(q))_ 

q , rn 

T 

L(q) = L(q,ia> m ,w„) 

= y ^ QG(Q+, iuJ m+ )G(Q-, iljJra). 

Q 

We see that for an isotropic system L(q)||q and the ver- 
tex corrections to the correlator of the orbital magnetic 
fields vanish exactly due to its particular vector struc- 
ture. 
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